A First-Principles Investigation of Phonon Softenings and Lattice Instabilities in the 

Shape-Memory System Ni2MnGa 
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Ferromagnetic Ni2MnGa has unique magnetoelastic properties. These are investigated by detailed 
computational studies of the phonon dispersion curves for the non-modulated cubic L2i and tetrag- 
onal structures. For the L2i structure, a complete softening of the transverse acoustic mode has 
been found around the wave vector q — [1/3, 1/3, 0](27r/a). The softening of this TA2 phonon mode 
leads to the premartensitic modulated super-structure observed experimentally. Further phonon 
anomalies, related to other structural transformations in Ni2MnGa, have also been found and ex- 
amined. These anomalies appear to be due to the coupling of particular acoustic phonon modes 
and optical modes derived from Ni. 



Displacive, diffusionless structural transformations of 
the martensitic type are known to occur in many metallic 
alloys. These transformations involve cooperative rather 
than diffusive displacements of atoms and are often as- 
sociated with phonon anomalies in the parent phase and 
related precursor phenomena. It is a general challenge of 
the fundamental physics to explain driving forces of the 
martensitic transformations. 

Since the discovery of a martensitic transformation in 
ferromagnetic NiaMnGa 0], this material has attracted 
strong interest jjl- This Heusler alloy is one of the 
very rare magnetic materials which undergo a martensitic 
transformation below the Curie temperature, whereby 
the combination of magnetic and structural features is 
responsible for its unique magnetomechanical properties. 
Shear deformations of more than 5 % have been ob- 
tained in magnetic fields Q These features to- 
gether make Ni2MnGa very efficient for magnetic-shape- 
memory (MSM) technology '6^. The MSM technology 
is based on the magnetic field induced redistribution of 
martensitic domains in the sample. From the techno- 
logical point of view, Ni2MnGa is much more promising 
than other materials being presently in commercial use, 
for example the well-known material Tb-Dy-Fe (Terfenol- 
D) which exhibits magnetostrictive strains of about 0.1 
%. Design of new efficient MSM magneto-mechanical ac- 
tuator devices is in progress • Also the search for new 
materials with magnetic shape-memory effects is under- 
way, even in antiferromagnets 0. 

Despite the experimental and technological success, a 
microscopic theory is missing, which would be able to 
show how far the premartensitic phase transition will re- 
sult from the coupling of soft phonon modes and homo- 
geneous strains associated with the shear constant C" 
I3, 01 • A complete understanding requires the evalua- 
tion of the phonon spectrum. The aim of this letter is to 
present supercell phonon calculations for the L2i and T 
structures (see Fig. QJand Tabled using a state-of-the- 



art first-principles method based on density-functional 
theory. 
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FIG. 1: (a) Simple cubic cell of Ni2MnGa (Heusler L2i struc- 
ture); (b) tetragonal bet cell in the [110] direction. The equi- 
librium lattice parameter obtained is flcubic = 5.8067 A and 
agrees with the experimental value Q. The direction [010]orth 
of the supercell corresponds to the direction [llOjcubic of the 
cubic structure. As the supercell consists of five tetragonal 
cells, the lattice parameters of the supercell are related to 
the lattice parameters of the cubic L2i structure as follows: 

Ctorth = (l/\/2)acubic, &orth = (5/\/2)acubic , Corth = Ocubic ■ 

Ni2MnGa is ferromagnetic at room temperature (Tc ~ 
380 K) and undergoes a two-step martensitic phase trans- 
formation. As a matter of fact, a number of different 
thermal and stress-induced martensitic structures have 
been observed in Ni2MnGa Tabled gives a sum- 

mary of the structures found for Ni2MnGa P, E, Oj- 
In the parent phase, several phonon anomalies can be 
observed with decreasing temperature. At the temper- 
ature T « 260 K, a nearly complete softening of the 
[1/3, 1/3, 0](27r/a) transverse acoustic TA2 mode with 
polarization along [110] leads to a premartensitic phase 
transition, which has been evidenced by neutron scat- 
tering |l3, 0) Q|> X-ray electron microscopy 
, and ultrasonic measurements [tR Il8l | or a combina- 



tion of the previous methods [l9j . This structural trans- 
formation involves a commensurate periodic distortion 
of the parent phase with a propagation vector equal to 
that of the soft mode q — [1/3, 1/3, 0](27r/a) correspond- 
ing to six atomic planes or three lattice spacings. The 
precursor phenomena involves the magneto-elastic cou- 
pling as it has been described using Landau-type models 
|2Q.] . These phenomenological models take as input very 
important lattice-dynamical properties. Thus, explicit 
calculation of the lattice dynamics is required since in 
transition metals, the d-electrons allow for many types 
of coupling, and several modes can be involved. 

TABLE I: Structural parameters of the crystal structures 
of Ni-Mn-Ga alloys according to Ref. [Tl| and this work. 
They follow (from left to right) the order of possible appear- 
ance when lowering the temperature. There is a superimposed 
modulation specified by the notation nL, where n stands for 
the close commensurate number of bet cells in the [110] di- 
rection; their respective q vectors are in units of {2n/a). The 
so-called premartensitic transformation happens between the 
L2i and 3L structures. 



Axes Lattice parameters (A) 





L2i 


3L 


5L 


7L 


T 


a 


5.824 


5.824 


5.90 


6.12 


6.44 


b 


5.824 


5.824 


5.90 


5.78 


5.52 


c 


5.824 


5.824 


5.54 


5.54 


5.52 


Modulation 


none 


q^ 0.33 


q ^ 0.43 


? 


none 


Tetragonafity 


1 


1 


0.94 




^ 1.2 



For the sake of simplicity, we present phonon calcula- 
tions for the direction [110], which is the most interest- 
ing one as seen from the measurements for the acoustic 
modes 0] . We have used the direct method for the eval- 
uation of the phonon dispersion curves [21I , whereby 
the Hellmann-Feynman forces are calculated with the Vi- 
enna Ab-initio Simulation Package (VASP) [ISIiil and 
the implemented projected augmented wave formalism 
(PAW) [i^. Within density- functional theory, the elec- 
tronic exchange and correlation are treated by using the 
generalized gradient approximation (GGA). The 3 c? elec- 
trons of Ga have been included as valence states. The 
importance of using PAW for Ni2MnGa instead of pseu- 
dopotentials has been pointed out earlier [2^. An or- 
thorhombic supercell of 40 atoms formed by five tetrag- 
onal crystallographic unit cells, as seen in Fig. ^b), has 
been used to calculate the phonon spectrum for the [110] 
direction jl^- The direct method implies that we cal- 
culate forces induced on all atoms of the supercell when 
a single atom is displaced from its equilibrium position. 
Displacement of the atoms in only one tetragonal unit cell 
along all directions allows one to derive the force constant 
matrix and the dynamical matrix. Diagonalization of the 
dynamical matrix then leads to a set of eigenvalues for 
the phonon frequencies and corresponding eigenvectors. 
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The calculations have been done by using an amplitude 
for the displacements oi u = 0.03 A. This value is sufB- 
cient to calculate the forces with required precision and 
to fulfill the conditions of the harmonic approximation. 

Phonon dispersion curves calculated for Ni2MnGa in 
the [110]cubic direction are shown in Fig. |2Ia). The 
acoustic modes are qualitatively and quantitatively in 
very good agreement with experimental results obtained 
from inelastic neutron scattering For instance, fol- 
lowing the acoustic branch, the values of the TA2 branch 
at C = 1, of the crossing between LA and TAi branch, 
and of the maximum of the LA branch are 2.33, 4.43 
and 5.13 THz, in good agreement with the experimental 
ones, 2.67, 4.35 and 5.07 THz (Ij. The initial slopes 
of the curves in the longitudinal {vl = 5.23 x 10^ m/s) 
and transverse (wtAi = 3.35 x 10^ m/s; VTA2 = 1.02 x 
10^ m/s) modes agree well with the sound velocities mea- 
sured via the neutron scattering dispersion curves |l3j . 
The slope of the phonon dispersion curves around the F 
point is positive, which is consistent with the stability 
around this phase. 




U110](27t/a) U110](27t/a) 



FIG. 2: Phonon dispersion curves for the cubic L2i structure 
(a) and the tetragonal T structure (b) of Ni2MnGa. Here, the 
reduced wave vector coordinate (" spans the fee BZ from F to 
X. 

We find a complete softening of the transverse acous- 
tic TA2 phonon mode between ^ = 0.2 and ( = 0.55. 
This softening occurs around the wave vector coordinate 
C — 1/3, which corresponds to the soft-mode phonon 
anomaly observed in the experimental studies [13. We 
point out that the experimentally observed premarten- 
sitic transition is related to this soft mode. 

The softening of the TA2 mode means that the cubic 
L2i high-temperature structure is unstable at zero tem- 
perature with respect to a particular atomic displacement 
leading to the formation of a modulated premartensitic 
structure. The 3L structure can be modeled by a large 
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FIG. 3: Polarization vectors corresponding to the acoustic 
phonon modes of Ga. The peaks at ^ = 0.300, ( = 0.433 and 
^ — 0.715 are due to the coupling to the optical modes of Ni. 



supercell composed of three basic bet supercells in the 
[110] direction (see Fig. Q 28,]. The atoms of the input 
structure have a modulated amplitude with a maximum 
around 0.05 A where the Ni atoms are in opposite phase 
to the Mn and Ga atoms. Full relaxation of the three 
lattice parameters, a, b and c has been allowed. Also, 
lower magnetizations of 60-80 % of the full magnetiza- 
tion value for the whole supercell has been taken into ac- 
count. In addition, two periods commensurate with the 
lattice were tested, but the structure with one period, as 
shown in Fig. 01 gives the minimum energy. The static 
displacements of the 3L modulated structure are shown 
in Fig. El^b). The modulation of the atoms stays in phase 
which differs from the input data, but it is typical for an 
acoustic branch such as TA2. The amplitudes for all the 
atoms are nearly the same, although they differ slightly 
but not according to the atom masses. Rather, the or- 
der is according to the sequence of the optical phonons 
energies, where Ni has the lowest energy, optical modes 
of Ga are in the next energy windows, and Mn optical 
phonons are of the highest energy (see Fig. [2J|. 

We have considered so far only the acoustic mode TA2 , 
because this mode is responsible for the phase transfor- 
mation L2i 3L. However, (see Fig. ^ the other two 
acoustic modes, TAi and LA, are coupled to the low - ly- 
ing optical modes of Ga. Fig. |3| shows how this coupling 
affects the behaviour of the polarization vectors of Ga for 
its acoustic modes. The Ga atom type have been chosen 
for clarity, while acoustic modes of the other atoms show 
similar trends. Figure O illustrates that the acoustic- 
optical hybridization (compare with Fig. [^J leads to a 
number of anomalies, which occur at the corresponding 
wave vector coordinates ( = 0.3, ( = 0.43 and ( = 0.715. 
A longitudinal-phonon optical mode can be associated to 
a charge density wave, which gives a sharp peak in the 
magnetic susceptibility x{q)- The crossing at C = 0.43 
compares well with those for which anomalies are ob- 
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FIG. 4: Displacements of the 3L premartensitic modulated 
structure: before a) and after b) structural relaxation. The 
atoms are displaced parallel to Oorth in successive (OlO)orth 
planes. The displacements are in a different scale in order to 
be observed. 



served in the susceptibility calculations, which have been 
associated with the lattice instabilities in the 5L and 7L 
structures 2jl, _2QJ . These crossings correlate in a natural 
way with the change of the C44 elastic constant, which 
is, as it turns out, not related to the TA2 branch through 
the premartensitic transformation. It comes as a novelty 
that the involvement of the optical modes is required in 
order to understand the phase transformations in these 
alloys. 

It is interesting to note that the anomaly in the 
TA2 branch has been found around the soft mode 
[1/3, 1/3, 0](27r/a) in our calculations, which is the wave 
vector associated with the phonon anomaly of the pre- 
martensitic transition [l^ . Also the phonon branches on 
the [110] direction for the T structure are given in Fig. 
2(b). In this latter case, the branches TA2 and TAi be- 
come similar and the softening of the TA2 branch, present 
in the cubic phase, disappears. This finding is consistent 
with the stability of the T structur e 1311 and adds insight 
into interpreting the experiments |l9l ]. The sequence of 
the optical modes is different in this case and they are 
not degenerate at the F and X points. The instability 
of the cubic structure at zero temperature is the clue to 
understand how the martensite phase may nucleate. In- 
stead of twinnings and stacking faults, the modulation of 
the martensitic tetragonal structure of Ni2MnGa corre- 
sponds to a smaller structural change easier to accommo- 
date or to derive from the parent cubic structure. With 
respect to the electronic structure, the lowering of the 
symmetry allows for a split ting of the density of states 
peak at the Fermi level [ifll, , thereby decreasing the 
valence band energy contribution to the total energy. 
Thus, structural modulation and lowering of the elec- 
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tronic energy explains the 5L structure of Ni2MnGa |32| . 

On the other hand, the recent susceptibUity calcula- 
tions 29, 3Q.J have shown that the peak in the gener- 
alized susceptibility at q = [(^, C, 0](27r/a) for C = 0.33 
corresponds to a magnetization of about 60%. This su- 
ceptibility peak moves toward higher values of C as the 
magnetization increases: full magnetization leads to a 
peak at C = 0.4. However our calculations show no need 
for the extra magnetization hypothesis, since all the in- 
teresting q vectors are already present. This casts serious 
doubts on the role of magnetization in these intermedi- 
ate or premartensitic transitions. Secondly, this finding 
reinforces the arguments that the elastic entropy term 
drives the transition 's^ and the spin spirals do not show 
anomalies involved in the shuffling of the L2i structure 
|33 |. In view of these results, the intermediate and the 
martensitic transition seem to be really independent, al- 
though, of course, the L2i structure is necessary in order 
to have the intermediate transformation at ( — 0.33. 

In summary, we have calculated the phonon disper- 
sion for two structures of Ni2MnGa: the cubic L2i and 
the tetragonal T at zero temperature. The cubic struc- 
ture is found to be unstable with respect to a specific 
re-arrangement (shuffling) of the atoms which leads to a 
modulated super-structure 3L. We have calculated sta- 
bility of the 3L structure with high accuracy and found 
that it reflects perfectly our phonon calculations. Fur- 
thermore, we have shown that the subtle competition of 
low-lying optical modes of Ni and acoustic modes ac- 
count for the experimentally observed wave vectors, for 
which anomalies have been observed (3L, 5L and 7L 
structures). These calculations have been done by using 
first-principles methods, and are in very good agreement 
with experimental findings. 
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